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It has been proposed that cosmic acceleration or inflation can be driven by replacing the Einstein- 
Hilbert action of general relativity with a function f(R) of the Ricci scalar R. Such f(R) gravity 
theories have been shown to be equivalent to scalar-tensor theories of gravity that are incompatible 
with Solar System tests of general relativity, as long as the scalar field propagates over Solar System 
scales. Specifically, the PPN parameter in the equivalent scalar-tensor theory is 7 = 1/2, which 
is far outside the range allowed by observations. In response to a flurry of papers that questioned 
the equivalence of f(R) theory to scalar-tensor theories, it was recently shown explicitly, without 
resorting to the scalar-tensor equivalence, that the vacuum field equations for 1/72 gravity around 
a spherically symmetric mass also yield 7 = 1/2. Here we generalize this analysis to f(R) gravity 
and enumerate the conditions that, when satisfied by the function f(R), lead to the prediction that 
7=1/2. 

PACS numbers: 04.50.+h ; 04.25.Nx 
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I. INTRODUCTION 

The evidence that the expansion of the Universe is cur- 
rently accelerating [l|, Q suggests that the Universe is 
dominated by dark energy with a large negative pressure. 
The predominant hypothesis is that a nonzero vacuum 
energy drives the acceleration, but this poses two serious 
theoretical questions: why is the vacuum energy nonzero, 
and why is it so miniscule? An equally plausible alterna- 
tive to dark energy is a modification of general relativity 
that would generate cosmic acceleration @,||. Modify- 
ing general relativity in this manner eliminates the need 
for dark energy, but it does not explain why the vacuum 
energy is zero. Similar modifications of general relativity 
have also been proposed to drive inflation || . 

A possible modification to general relativity that gen- 
erates an accelerated expansion is 1/i? gravity ||, in 
which a term proportional to 1/i?, where i? is the Ricci 
scalar, is added to the Einstcin-Hilbert action so that 
the 1/i? term dominates as the Hubble parameter de- 
creases. Soon after the introduction of this theory, it 
was shown that 1/i? gravity is dynamically equivalent 
to a scalar-tensor gravity with no scalar kinetic term 
Q. Moreover, the equivalence to scalar-tensor gravity 
applies to all modified gravity theories that replace the 
Einstcin-Hilbert action with some function of the Ricci 
scalar [known as f(R) gravity], provided that f(R) has a 
nonzero second derivative with respect to R. When the 
scalar field is light, this theory makes predictions that are 
incompatible with Solar System tests of general relativity 
@, H, 9 ■ Consequently, Ref . [|| concluded that a broad 
class of /(i?) gravity theories, including 1/i? gravity, are 
ruled out by Solar System tests. 

Since then, however, the results in Ref. 
cized by a number of papers [lj| EH El, EH 



even claim that Solar System experiments do not rule 
out any form of /(i?) gravity. The essence of the crit- 
icism is that /(i?) gravity admits the Schwarzschild-dc 
Sitter solution and hence the vacuum spacetime in the 
Solar System is not different from that in general rela- 
tivity, although there were also broader objections to the 
equivalence between /(i?) and scalar-tensor gravity [u|. 
Working directly with the field equations, a recent paper 
[IH found that even though the Schwarzschild-de Sitter 
metric is a vacuum solution in 1/i? gravity, it does not 
correspond to the solution around a spherically symmet- 
ric massive body. 1 They found that the solution for the 
Solar System is identical to the spacetime derived using 
the corresponding scalar-tensor theory 

In this paper, we generalize the analysis of Ref. [l5j to a 
broad class of f(R) gravities, namely those theories that 
admit a Taylor expansion of /(i?) around the background 
value of the Ricci scalar. We work in the metric formal- 
ism, where the field equations are obtained by varying 
the action with respect to the metric and treating the 
Ricci scalar as a function of the metric. The Palatini 
formalism, which treats the Ricci scalar as a function of 
the connection and varies the action with respect to the 
connection and the metric independently, yields differ- 
ent field equations for /(i?) gravity and has been studied 
extensively elsewhere (e.g. Refs. 0, E US HI H ) ■ 

This paper is organized as follows: In Section HH we 
solve the linearized field equations around a spherical 
mass and find that the solution in the Solar System is 
in agreement with the solution obtained using the equiv- 



were criti- 
I and some 



1 Eddington made a similar mistake in R 2 gravity lid , which 
later corrected by |17| . 
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alent scalar-tensor theory. When f(R) satisfies a condi- 
tion that is analogous to the scalar field being light in the 
equivalent scalar-tensor theory, the resulting spacetime is 
incompatible with Solar System tests of general relativ- 
ity. In Section |III[ we consider how our analysis applies 
to several f(R) gravity theories, including general rela- 
tivity. This particular example illustrates the connection 
between f(R) gravity and general relativity and clarifies 
the requirements for a general relativistic limit of an f(R) 
theory. We summarize our conclusions in Section IIVI 



II. 



two components: 



WEAK-FIELD SOLUTION AROUND A 
SPHERICAL STAR 



We consider gravitational theories with actions of the 
form 



1 

2^7 



(1) 



where f(R) is a function of the Ricci scalar R and S m is 
the matter action. The field equation obtained by vary- 
ing the action with respect to the metric is 

IrR^u - \^!g^v - V^Vjz/r + OfRg^u = KTp,,, (2) 

where f R = df/dR. In previous studies, predictions of 
Solar System dynamics in these theories were analyzed 
by appealing to an equivalence with scalar-tensor theo- 
ries |6(. We review this equivalence in Appendix^ Since 
the equivalent scalar-tensor theory is incompatible with 
Solar System observations if the scalar field propagates 
on Solar System scales, Ref. @ concluded that the cor- 
responding f(R) theories are ruled out. We now show 
that this conclusion can be made without appealing to 
the equivalence between f(R) and scalar-tensor gravity. 
Instead, we work directly with the linearized field equa- 
tions about a spherical mass distribution. Our treatment 
clarifies and amends a similar analysis presented in Ref. 
p3 |. and we extend it to cases where the background 
value of the Ricci scalar equals zero. 

We now find the metric that describes the spacetime 
around a spherical body in f(R) gravity in the weak-field 
regime. To do this, we must choose a background space- 
time around which to linearize the field equations. The 
only physically relevant choice is an isotropic and ho- 
mogeneous background spacetime that solves Eq. @ for 
some spatially uniform cosmological stress-energy tensor 
T^° s . The evolution of the time-dependent and spatially 
homogeneous background scalar curvature Ro(t) is deter- 
mined by the trace of Eq. ©, 

fRo(t)Ro(t) - 2/o(«) + 3D/j«,(t) = KT COB (t), (3) 
where f m = df/dR\ R=Ra , f = f(R Q ) and T cos = 

y ■ 

In order to investigate perturbations away from this 
background, we express the Ricci scalar as the sum of 



R(r,t) = R (t) + Rxir), 



(4) 



where Ro(t) is the spatially homogenous background 
curvature that solves Eq. ([3]) and Ri(r) is a time- 
independent perturbation to this background curvature. 
We assume that all derivatives of f(R) arc well-defined 
at the present-day value of Rq so that we may use a 
Taylor expansion of f(R) around R = Rq to evaluate 
/(i?o + Ri) and f R (Ro + Ri). We will terminate the ex- 
pansion by neglecting terms nonlinear in R\. Provided 
that the higher-order terms of the Taylor series do not 
cancel in some contrived way, neglecting the higher-order 
terms is only justified if the sum of the zeroth-order and 
linear terms is greater than all other terms in the Taylor 
expansion. Specifically, we require that 

f0 + fR0Rl » V 1 "'*)^' ( 5 ) 

n! 

fRO + fRRoRi » ^f {n+1) (Ro)Ri, for all n> 1,(6) 
n! 

where f RR0 = d 2 f/dR 2 \ R=Ro and f^(R ) = 
d n f/dR n \ R = Ro . 

Now we consider the trace of Eq. ([2|) with both a 
cosmological matter source described by T cos and a finite, 
time- independent, spherically symmetric matter source, 
described by T s : 



f R R - 2/ + 3D/ fl = K (T cos + T s ) . 



(7) 



Using first-order Taylor expansions to evaluate f R and 
/ and neglecting 0(R\) terms, we obtain a linearized 
version of Eq. ([7]) : 



3/ M0 n-Ri(r) 



-fRR0(t)R (t) ~ 3D/flflo(t) 



i?i = kT s 



(8) 



To obtain this equation, we used the fact that Ro(t) 
solves Eq. ([3]) to eliminate terms that are independent 
of R\. By dropping 0(f RR oR 2 ) terms from Eq. while 
keeping the f RR oRoRi term, we have implicitly assumed 
that i?i <C i?o if Ro is nonzero. We will check that 
this condition is satisfied after the discussion following 
Eq. (|3"Tj) . If R is zero, then the 0(f RRa R\) is guaran- 
teed to be smaller than the nonzero terms in Eq. (|5J) by 
virtue of Eq. ([5]). Note that if f RR o = 0, as in general 
relativity, this equation becomes simply f R aRi = —kT s . 
If in addition f R o is nonzero then R\ must vanish outside 
the star and hence the Schwarzschild-de Sitter solution 
becomes the solution to the field equation outside the 
source. However, if f RR o ^ 0, this is no longer necessar- 
ily the case. 

Finally, we take our background metric to be a flat 
Fricdmann-Robertson- Walker (FRW) metric. We then 
consider a spherically symmetric perturbation to this 
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background so that the linearized perturbed metric takes 
the form 



ds 2 



[l + 2f(r)]dt 2 
-a(tf{ [1 + 2$(r)]dr 2 + r 2 aV? } , 



(9) 



where the present value of a(t) is one. When solving 
the field equations, we will keep only terms linear in the 
perturbations \& and 

We will now solve Eq. (jSJ) for a nonzero fnm- Since 
we confine our analysis to a static perturbation R\{r), □ 
becomes the flat-space Laplacian operator V 2 . Restrict- 
ing our analysis to a source with mass density p(r) and 
negligible pressure, we may rewrite Eq. (JSJ) as 



V 2 i?i - m 2 R 1 



up 



where we have defined a mass parameter 



1 ( I- 



RO 



3 \fllR0 



Ro 



Pf- 



rro 



Jrbo 



(10) 



(11) 



Due to the evolution of Rq (t) , this mass parameter varies 
in time. However, the time-scale of variation in the cos- 
mological background spacetime is comparable to the 
current Hubble time. Since this time-scale is much longer 
than the time-scale of Solar System dynamics, we may 
neglect the time variation of the background spacetime 
when considering the behavior of bodies within the Solar 
System [24[. Therefore, for the purposes of this calcula- 
tion, we take m to be time-independent. 

The Green's function G(r) for this differential equation 
depends on the sign of m 2 : 



G(r) 



— cos(mr) / (Anr) m 2 < 0, 

— exp(— mr) / (4irr) m 2 > 0, 



(12) 



where m = y/\m 2 \. If mr <C 1, then both Green's func- 
tions are approximately — l/(47rr), which is the Green's 
function for Laplace's equation. In this case, the term 
proportional to m 2 in Eq. (fTU|) may be neglected and the 
solution outside the star is given by 



Ri = 



M 



12nf RR0 r 



(13) 



where M is the total mass of the source. We note that 
when applied to l/R gravity with a static de Sitter back- 
ground, this result agrees with the result presented in 
Ref. 

We emphasize that in order for this solution for Ri 
to be valid, we must have mr <C 1. Only when this 
condition is satisfied is the trace of the field equation 
well-approximated by Laplace's equation. This restric- 
tion was not mentioned in Ref. [23[ ■ The physical inter- 
pretation of this constraint is clear when one considers 
the equivalent scalar-tensor theory. When one switches 
to a frame where the scalar degree of freedom is canoni- 
cal, the effective mass of the scalar field evaluated in the 



Jordan frame is | 
2 fao ( 



RRO 



Rq 
/ro 



4/o 



2kT c 



(Jro) 2 (Jro) 2 



(14) 



Since Rq is the solution to Eq. ([3]), this expression may 
be simplified to 



RO 



-Jin -6 



/,'(> 



RRO 



RO 



(15) 



It is clear that both m v and m [defined by Eq. (fTTj) ] are 
of the same order. Therefore, the condition that mr <C 1 
is equivalent to demanding that the scalar field be light 
(m v r -C 1). See Appendix A for more details. 



In summary, Eq. (|T3|) is a solution to the trace of the 
field equation within the Solar System only if the scalar 
degree of freedom propagates on Solar System scales. In 
terms of f{R), the necessary condition is 



ImV = 



RO 



— Rn — 3 



RRO 



RRO 



RRO 



r 2 < 1. 



(16) 



The triangle inequality tells us that the mass constraint 
given by Eq. (fl6|) implies that 



RO 



fRRO 



Ro 



PI 



RRO 



fRR.O 



< 1. 



(17) 



Finally, since Of RR0 / f Rm ~ H 2 , where H = a/a is the 
current Hubble parameter, and we know that Ror 2 ~ 
H 2 r <C 1 by cosmological constraints, the mass con- 
straint implies that 



/fto 



RRO 



< 1. 



(18) 



We will now use the expression for R\ given by Eq. (JTJ 
to solve the field equations for the metric perturbations 
^ and <f>. As we did for the trace of the field equa- 
tion, we simplify the field equations by replacing f(R) 
and f R (R) with first-order Taylor expansions around the 
background value Rq to obtain field equations that are 
linear in R\. Using Eq. ([3]) to simplify this expression, 
we obtain 



fm{Ri — \Rq\v) + JrrqRiRv ~ -j^fRoRiK 
-fRnoVVvRi + SSfiuuClRi 



(19) 



where [Ro] 1 / is the unperturbed FRW Ricci tensor and 
5% is the Kronecker delta. We neglected time derivatives 
of the background metric when deriving this equation. 
As previously noted, the time-scale of variations in Rq is 
much longer than that of Solar System dynamics, making 
the terms involving time derivatives of Rq irrelevant to 
gravitational effects within the Solar System. 

We simplify Eq. (|19[) further by dropping several negli- 
gible terms. We continue to ignore terms that depend on 
the variation of the background spacetime by dropping 
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terms that involve products of $, \P and JhroRi with 
H and dH/dt. Since we are working in the weak-field 
regime, we neglect all terms that are nonlinear functions 
of the metric perturbations $ and \&. Keeping only terms 
that are linear in $ and iff allows us to replace the □ 
with the flat-space Laplacian operator V 2 since the per- 
turbation is assumed to be static. Finally, we know from 
Eq. (JT3J) that /rroRi ~ kM/t, and we expect \I> and <E> 
to be proportional to kMJt as well. Therefore, /rro-Ri^ 
and fRRoRi® are second-order quantities, and we may 
neglect them. With these simplifications, the tt,rr,96 
components of Eq. (TIT)]) are respectively 

fm (-*" + - IfnoRi + -fRRvR'i = 0, (21) 
\ r J 2 r 

{ 1 1 2 

\r r r z 

— T^fnoRi + -fRRoR'i + IrroR'i = 0, (22) 

where the prime denotes differentiation with respect to 
r. The (fxp component of Eq. I|19p is identical to the 09 
component given by Eq. (|22|) . 

Recalling that R\ solves Eq. (fT0|) with m 2 = so that 
V 2 i?i is proportional to the density p, Eq. ([20]) may be 
rewritten 



1 



f R oV^=-Kp--f RO Ri 



(23) 



We express iff as the sum of two functions: iff = Wo + 
where 



/floV 2 *! = -i/Roiii. 



(24) 
(25) 



Provided that f R o ^ 0, Eq. (|24p may be integrated via 
Gauss's Law to give 



k m(r) 
671-/^0 r 2 



(26) 



where m(r) is the mass enclosed in a sphere of radius 
r. If we assume that "to vanishes as r — > oo, we may 
integrate Eq. (|26|) to obtain 



71/ 



(27) 



outside the star. Solving Eq. (J25J) outside the star using 
Eq. (H3D for i?! yields 



1 



^8-kJrro 



nMr < 



1 kM 
fm r 



(28) 



where the inequality follows from Eq. (jTHJ) • Since Wo ~ 
kMKJrqt) outside the star we have shown that |Wi| <C 



|Wo|. Therefore, we may neglect Wi and conclude that 
f = $o ^ given by Eq. P7|) . This expression for iff is 
used to define Newton's constant: G = K/{6nfm). For 
1/R gravity with a static vacuum de Sitter background, 
fm = 4/3, so k takes its standard value of 8irG and 
Eq. I|27p matches the corresponding result in Ref. [l5[. 

We now turn our attention to Eq. (|2ip . which we will 
solve for <!>. First, we note that Eq. (| 13(1 implies that 
R[ = —Rx/r. Therefore, the ratio of the second two 
terms in Eq. (|2~lj) is 



(l/2)/floi?i 




JflO 


VRmR'Jr 




fRRO 



r z < 1, 



(29) 



where the inequality follows from Eq. (|18p . Consequently, 
the fRoRi term is negligible, and we drop it from the 
equation. Differentiating Eq. (f26|) to find iff", and using 
Gauss's Law to obtain R[ from Eq. pTT)) (with m 2 = 0), 
we may then rewrite Eq. (|2"Tj) as 



$'(r) 



k d / m(r) 



12nf R0 dr 



(30) 



Assuming that $ vanishes as r — > oo, this equation may 
be integrated to obtain 



A/ 



127r/ fl0 r ' 



(31) 



outside the star. It is easy to verify that Eqs. (|2?[) and 
([31) also satisfy the third field equation, Eq. ([22]) . 

We may now check our assumption that R\ -C i?o for 
nonzero i?o- From the expression for R\ given by Eq. (j!3|) 
and our definition that k = GirfRoG, we see that 



0±£ — ( GM \ f«o 
Ro Rq \ R s J f_ 



RRO 



(32) 



where R s is the radius of the star. It is easy to check that 
this expression holds inside the star as well by integrating 
Eq. ([TU[) into the interior of the star. Therefore, our 
assumption that R\ -C Ro places an additional condition 
on the ratio fR /fitm- 



f 



R0 



fRRO 



< Ro 



\GM J 



for Rq 0. 



(33) 



If fRo/ fRRO ~ Ro, as is the case for many f{R) theories 
with nonzero Ro, then this condition is always satisfied. 

Thus we have shown explicitly that iff = —2$ = 
—GM /r for all f(R) theories with nonzero fRRo that sat- 
isfy the conditions given by Eqs. (O, ([6]), (fT6[) and ([33]) . 
Transforming the metric given by Eq. ([9]) to isotropic co- 
ordinates, taking a — 1 today, and keeping only terms 
that are linear in GM/r gives 



ds 2 



1 - 



+ 1 



2GM\ , 
—) dt 



GM 



[dr 



'-dn 2 



(34) 
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It is clear that this spacctimc is equivalent to a Parame- 
terized Post-Newtonian spacetime with PPN parameter 
7 = 1/2. This result is in gross violation of observations; 
Solar System tests require that 7 = 1 + (2.1 ± 2.3) x 10~ 5 
We also note that this result is in precise agreement 
with the results obtained using the equivalent scalar- 
tensor theory [|| (see also (25j). 



III. CASE STUDIES 

First, we show how we regain the results of general 
relativity if we take fnno — and assume that our lin- 
earized Taylor expansion is a valid approximation. We 
note that general relativity [f(R) = R] satisfies both of 
these conditions. 

Taking f RR0 = 0, Eq. © yields 



frtoRi — *>P- 



(35) 



When Jrro = 0, the JrqRi terms in the field equations 
[Eqs. (|21ti22|) ] are no longer negligible compared to the 
terms proportional to Jrro since these terms vanish. The 
field equations then become 

/aoV 2 * + i/ fl ofli = k P , (36) 

fno + - \faoRi = 0, (37) 

/*> ( - + - h R oRi = 0. (38) 



Using Eq. ((35]), Eq. ([36]) becomes 



fmV 2 ^ = -p, 



and the solution outside the star is 

K M 

* = 



87r//{o r ' 
From Eq. (J37J) and Eq. ([38]), we have 



and the solution outside the star is 
« M 



(39) 



(40) 



(41) 



(42) 



Since * = — $ = —GM/r, transforming to isotropic co- 
ordinates reveals that 7 = 1 as expected. 

With this result it is easy to see why the p — > limit in 
1 /R n (n > 0) gravity does not recover general relativity. 
In 1/R n gravity [|[, we have 



The static solution to Eq. © with T cos = is R = (n + 
2) 1 /(" +1) /i 2 , and f RR0 cx ^ 2 . Therefore, f RR0 diverges 
rather than vanishes in the limit that fx — > 0, and general 
relativity is not regained. The mass parameter for this 
theory has the dependence m 2 cx jj 2 and hence it vanishes 
in the limit that p, — ► 0. Furthermore, a Taylor series 
of Eq. (|43|) around i?o is well-behaved and cosmological 
constraints tell us that /i ~ so that m 2 r 2 <C 1 in the 
Solar System. We conclude that the analysis of general 
f(R) gravity given in Section UT1 applies and 7 = 1/2 for 
these theories in a static background. 

We note however that the static solution to Eq. ([3]) 
may not describe the current cosmological background 
in l/R n gravity. This solution is unstable, and without 
fine-tuning of the initial conditions, this spacetime will 
evolve toward a spacetime with Rq <C p 2 0|. In that 
case, we note that 



(m\)- 1 f ( - m \R )R^ 
jo + JroRi 
(m\)- 1 f < - m +^{R )R^ 
fm + JrroRi 



< 



< 1, 



« 1, 



(44) 
(45) 



so that Eqs. ([5]) and ([6]) are still satisfied. Furthermore, 
m 2 cx i?o, so, as in the static-background case, the mass 
is of order the Hubble parameter today. Therefore, the 
7=1/2 result holds even during the late-time evolution 
of 1/R n gravity. 

Next we consider Starobinsky gravity Q which has 



(46) 



The static solution to Eq. © with T cos = is R Q = 
for this theory. Since f(R) is a second-order polynomial, 
the first-order Taylor expansion of /i?(i?o + Ri) is exact. 
The 0{R 2 ) term in the Taylor expansion of f(Ro + Ri) 
is suppressed compared to the linear term by a factor of 
GM/r and is therefore negligible. The mass parameter 
for this theory is proportional to a 2 , so Eq. (TIB"]) is a 
solution for i?i if a 2 r 2 -c 1. Therefore, 7 = 1/2 in 
this theory if a 2 r 2 <C 1 inside the Solar System. If the 
mass parameter a is made large (i.e. if a ~ 10 12 GeV as 
proposed in Ref. [f|), then this condition is not satisfied 
and we cannot use the analysis in Section [IT] to calculate 
7 for this theory. 

Next we consider an example of a theory that uses two 
mass parameters: a hybrid between Starobinsky gravity 
and l/R gravity. In particular, consider the function 



f{R) = R+- 2 R 2 -^- 
a z R 



(47) 



We then find that, as in the usual l/R case, we have i?o 
\J~Sp 2 (for a static background in vacuum). However, 



f(R) = R 



2+2)1 



R" 



n > 0. 



(43) 



2 o 2 

m = 3/x 



(48) 
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We can make this quantity as large as we want by let- 
ting the denominator tend towards zero, which gives the 
condition a — > 3 3//4 /i. Thus, in this model we can vio- 
late the conditions listed in Section [IT] by fine-tuning the 
parameters. 

Finally, we consider power-law gravitational actions 



f(R) 



1+5 



(49) 



Assuming that 8^1, the static vacuum solution to 
Eq. is Rq = 0. If 5 is not an integer, there will 
be some derivative that is not defined at R = 0, which 
causes the Taylor expansion to fail around that point. In 
particular, if it is supposed that 5 <^ 1, then at least the 
second derivative will be undefined so that the Taylor ex- 
pansion will fail. For 5 = 1 the static vacuum background 
value i?o is undetermined. However, if we choose Rq ^ 
then all of the conditions listed in Section [IT] are satis- 
fied and we conclude that 7 = 1/2 in agreement with 
Ref. [27[. If 8 is an integer greater than one, then the 
Taylor expansion around f(Ro = 0) is well-defined, but 
we cannot drop the terms that are nonlinear in R\ since 
the linearized function vanishes. Therefore, this analysis 
is incapable of determining whether f(R) = R 1+5 gravity 
with 5^1 conflicts with Solar System tests. 



IV. CONCLUSIONS 

By analyzing the field equations around a spherically 
symmetric mass, we have shown that, in agreement with 
the analysis in Ref. @ , the PPN parameter 7 of general 
f(R) gravity is 7 = 1/2 given the following conditions: 

I. The Taylor expansions of f(R) and df/dR about 
the current background value R = Rq, where Rq solves 
Eq. ([3]), are well-defined and dominated by terms that 
are linear in deviations away from R = Rq. If Rq is non- 
zero, then the deviations from Rq are small compared 
to i?o- This condition may be re-expressed as Eq. (j33|) 
and is closely related to the third condition stated below. 

II. The second derivative of f(R) with respect to 
R is nonzero when evaluated at the background value of 
R = Rq. 

III. The mass parameter given by Eq. (fTTj) respects the 
condition mr <C 1 within the Solar System. 



For theories with one extra mass parameter and 
non-zero Rq, as in l/R gravity, it is reasonable to 
assume that fao/faRo ~ Rq. In that case, the latter 
part of the first condition is always true and the third 
condition is satisfied provided that Rqt 2 -C 1 within 
the Solar System. However, for theories with multiple 
mass parameters, such as the Starobinsky-l/i? hybrid 
presented in this paper, it is possible that this condition 
can be violated. 

The second and third conditions listed above corre- 
spond to synonymous conditions in the scalar-tensor 
treatment: f(R) and scalar-tensor gravity are equiva- 
lent only if the second derivative of f(R) is nonzero, 
and 7 = 1/2 only if the scalar field is light enough to 
propagate through the Solar System. Therefore, we have 
also verified that, contrary to the claim of some authors 
fiol . [Til IT2I IT3L [Tj] , calculating the Solar System predic- 
tions of f(R) gravity using the equivalent scalar-tensor 
theory is a valid technique. 
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APPENDIX A: REVIEW OF SCALAR TENSOR 
EQUIVALENCE 

The action for the scalar-tensor theory that is equiva- 
lent to f(R) gravity is 



S 



1 

2k 



dW^[/(0) + / (0)OR-0)] + S m , (Al) 



where /</,(</>) = df /dip and S m is the matter action. The 
field equation for <fi is <j> = R if d 2 f/d<fr 2 ^ 0. Since 
the relation between </> and R is purely algebraic, it can 
be resubstitutcd into the action to reproduce the action 
for f(R) gravity given by Eq. |T]). After the conformal 
transformation = f4,(<f>)g^, the action becomes that 
of general relativity with a minimally coupled scalar field: 



S = 7^7 J d i x y /-g E (^Re 



r<?nV S/ J^)][VW0(0)] 



[MM - m\ + z 



(A2) 



7 



Introducing a canonical scalar field ip such that f$(4>) = 
cxp(-\/2k/3</j), Eq. (|A2[) can be rewritten as 

S = J cfxV^^RE-^iVEV) 2 ~V(ip)j +S m , 

(A3) 

where the potential is defined by 



V{<p) 



(A4) 



The absence of the kinetic term in Eq. (|A1|) implies the 
Brans-Dickc parameter of f(R) gravity theories is u> = 
■ From an analysis of Brans-Dicke gravity, if the scalar 
degree of freedom can propagate on scales much larger 
than the Solar System, we can conclude that 7 = (1 + 
w)/(2 + w) = 1/2 Q. 

In the frame where p> is canonical (the Einstein frame) 
ip has the equation of motion 



(A5) 



where the prime denotes differentiation with respect to 
4>. When we re-express Eq. (| A5|) in terms of f'(4>) and 



the usual metric g fiVl we recover Eq. ((7]). Therefore, we 
stress that this reformulation contains no new dynamics 
compared to the expressions used in this paper. The two 
formulations are entirely equivalent. 

In order to derive the mass m v , we let p = <pa(t)+ipi (r) 
and T M = T cos + T s so that ip (t) satisfies Eq. CM} with 
T cos . We then expand to linear order in the perturbation 
tpi, writing Eq. (|A5|) in terms of the physical metric g^ v . 
We find 



□pi = f'fo) 



d 2 V 



dp 2 



3 [/'(</>o)] 5 



K T s 



6 f'(M ' 
(A6) 

where 4>o denotes the background field value for the <fi 
field. Using Eq. (|A4j) to evaluate d 2 V/dtp 2 , we have 



1 



00 4/((/) ) 



2k 



into)) 2 . 

(A7) 



Finally, we may rewrite m 2 as Eq. ()14[) since tfio = Rq. 
We conclude that if m^r 2 <C 1 then 7 = 1/2 as discussed 
in Ref. 0. 
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